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ABSTRACT 



Ongoing searches for supernovae (SNe) at cosmological distances have recently started 
to provide a link between SN la statistics and galaxy evolution. We use recent estimates 
of the global history of star formation to compute the theoretical Type la and Type 
II SN rates as a function of cosmic time from the present epoch to high redshifts. We 
show that accurate measurements of the frequency of SN events in the range < z < 1 
will be valuable probes of the nature of Type la progenitors and the evolution of the 
stellar birthrate in the universe. The Next Generation Space Telescope should detect of 
order 20 Type II SNe per 4' x 4' field per year in the interval 1 < z < 4. 



Subject headings: galaxies: evolution - supernovae: general 
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1. Introduction 

The remarkable progress in our understanding of faint galaxy data made possible by the 
combination of HST deep imaging (Williams et al. 1996) and ground-based spectroscopy (Lilly et 
al. 1995; Ellis et al. 1996; Cowie et al. 1996; Steidel et al. 1996), has recently permitted to shed some 
light on the evolution of the stellar birthrate in the universe, to identify the epoch 1 < z < 2 where 
most of the optical extragalactic background light was produced, and to set important contraints 
on galaxy formation scenarios (Madau et al. 1998; Steidel et al. 1998). While one of the biggest 
uncertainties in our knowledge of the emission history of the universe is probably represented by 
the poorly constrained amount of starlight that was absorbed by dust and reradiated in the IR 
at early and late epochs, one could also imagine the existence of a large population of relatively 
old or faint galaxies still undetected at high-z, as the color-selected ground-based and Hubble Deep 
Field samples include only the most actively star-forming young objects. It is then important 
at this stage to devise different observational strategies, free of some of the biases that plague 
current galaxy surveys, and to make testable predictions for future astronomical capabilities, such 
as SIRTF, FIRST and NGST. 

Here, we shall focus our attention on the rate of supernova (SN) explosions in the universe. 
An obvious reason to consider SNe is purely observational, i.e. the fact that they are very bright 
objects, with luminosities as high as 10 10 L@, and are point-like sources, making their detection 
possible even at very large distances/redshifts. More in general, the evolution of the SN rate with 
redshift contains unique information on the star formation history of the universe, the initial mass 
function (IMF) of stars, and the nature of the binary companion in Type la events. All are essential 
ingredients for understanding galaxy formation, cosmic chemical evolution, and the mechanisms 
which determined the efficiency of the conversion of gas into stars in galaxies at various epochs 
(e.g. Madau et al. 1996; Madau, Pozzetti, & Dickinson 1997; Renzini 1997). While the frequency 
of "core-collapse supernovae", SN II and possibly SN Ib/c, which have short-lived progenitors 
(e.g. Wheeler & Swartz 1993) is essentially related, for a given IMF, to the instantaneous stellar 
birthrate of massive stars, Type la SNe - which are believed to result from the thermonuclear 
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disruption of C-0 white dwarfs in binary systems - follow a slower evolutionary clock, and can 
then be used as a probe of the past history of star formation in galaxies (e.g. Branch et al. 1995; 
Ruiz-Lapuente, Canal, & Burkert 1997; Yungelson & Livio 1998). The recent detection of Type 
la SNe at cosmological distances (Kim et al. 1997; Garnavich et al. 1998; Perlmutter et al. 1998) 
allow for the first time a detailed comparison between the SN rates self-consistently predicted by 
stellar evolution models that reproduce the optical spectrophotometric properties of field galaxies, 
and the observed values. 

In this Letter we show how accurate measurements at low and intermediate redshifts of the 
frequencies of Type II(+Ib/c) and la SNe could be used as an independent test for the star formation 
and heavy element enrichment history of the universe, and significantly improve our understanding 
of the intrinsic nature and age of the populations involved in the SN explosions. A determination 
of the amount of star formation at early epochs is of crucial importance, as the two competing 
scenarios for galaxy formation, monolithic collapse - where spheroidal systems formed early and 
rapidly, experiencing a bright starburst phase at high-z (Eggen, Lynden-Bell, & Sandage 1962; 
Tinsley &: Gunn 1976) - and hierarchical clustering - where ellipticals form continuosly by the 
merger of disk/bulge systems (White &: Frenk 1991; Kauffmann et al. 1993) and most galaxies 
never experience star formation rates in excess of a few solar masses per year (Baugh et al. 1998) - 
appear to make rather different predictions in this regard. We show how, by detecting Type II SNe 
at high-z, the Next Generation Space Telescope should provide an important test for distinguishing 
between different scenarios of galaxy formation. 

2. Basic Theory 

2.1. Cosmic Star Formation History 

In this section we shall follow Madau et al. (1998), and model the emission history of field 
galaxies at ultraviolet, optical, and near-infrared wavelengths by tracing the evolution with cosmic 
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time of their luminosity density, 



DO 




dL u L v y{L v , z) = T(2 + a)y*L 



(1) 



Jo 



where y(L u ,z) is the best-fit Schechter luminosity function in each redshift bin. The integrated 
light radiated per unit volume from the entire galaxy population is an average over cosmic time 
of the stochastic, possibly short-lived star formation episodes of individual galaxies, and follows a 
relatively simple dependence on redshift. Madau et al. (1998) have shown how a stellar evolution 
model, defined by a time-dependent star formation rate per unit volume, ip(t), a universal IMF, 
4>{m), and some amount of reddening, can actually reproduce the optical data reasonably well. In 
such a system, the luminosity density at time t is given by the convolution integral 



where l v (t') is the specific luminosity radiated per unit initial mass by a generation of stars with 
age t', and p esc is a time-independent term equal to the fraction of emitted photons which are not 
absorbed by dust, p esc = exp(— r u ) in a foreground screen model. The function ip{t) is derived from 
the observed UV luminosity density, and is then used as input to the population synthesis code of 
Bruzual & Chariot (1998). 

Figure la shows the model predictions for the evolution of p u for a Salpeter (1955) function 
(4> oc m -2 ' 35 ), .Eb-v = 0.1 with SMC-type dust (in this case, the observed UV luminosities must 
be corrected upwards by a factor of 1.4 at 2800 A and 2.1 at 1500 A), and a star formation history 
which traces the rise, peak, and sharp drop of the UV emissivity. ^ For simplicity, the metallicity 

2 A convenient fit to 5% accuracy to the star formation history used for Figure la is ip(t) = 
0.049[^ e - i9 /°- 64 + 0.2(1 - e-* 9 /°- 64 )]M yr- 1 Mpc~ 3 , where t 9 is the Hubble time in Gyr, t 9 = 
13(1 + z) -3 / 2 . Note that, although in our calculations the IMF extends down to mi = 0.1 Mq, 
stars below 0.8 Mq make only a small contribution to the emitted near-IR light. This introduces a 
non-negligible uncertainty in our estimates of the total stellar birthrate. For example, in the case of 
a Salpeter IMF with mi = 0.5 Mq, the inferred star formation rate, ip oc (m ; ~ ' 35 — m~ ' 35 ), would 
decrease by a factor of 1.88. 




(2) 
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was fixed to solar values and the IMF truncated at 0.1 and 125 M : none of these assumptions is 
crucial for our results. The data points show the observed luminosity density in six broad passbands 
centered around 0.15, 0.20, 0.28, 0.44, 1.0, and 2.2 fj>m. The model is able to account for the entire 
background light recorded in the galaxy counts down to the very faint magnitude levels probed by 
the Hubble Deep Field (HDF), and produces visible mass-to- light ratios at the present epoch which 
are consistent with the values observed in nearby galaxies of various morphological types. The bulk 
(<; 60% by mass) of the stars present today formed relatively recently (z < 1.5), consistently with 
the expectations from a broad class of hierarchical clustering cosmologies (Baugh et al. 1998), and 
in good agreement with the low level of metal enrichment observed at high redshifts in the damped 
Lyman-a systems (Pettini et al. 1997). 

One of the biggest uncertainties in our understanding of the evolution of luminous matter in 
the universe is represented by the poorly constrained amount of starlight that was absorbed by 
dust and reradiated in the far-IR at early epochs. While the Salpeter IMF, -Eb-v = 0.1 SMC- 
dust model reproduces quite well the rest-frame UV colors of high-z objects in the HDF, the 
prescription for a "correct" de-reddening of Lyman-break galaxies is the subject of an ongoing 
debate (see Pettini et al. 1998 and references therein). Figure 16 shows the model predictions 
for a monolithic collapse scenario, where half of the present-day stars - the fraction contained in 
spheroidal systems (Schechter & Dressier 1987) - were formed at z > 2.5 and were enshrouded by 
dust.[] Consistency with the HDF "dropout" analysis has been obtained assuming a dust extinction 
which increases rapidly with redshift, Eb-v = 0.011(1 + z) 2 ' 2 . This results in a correction to the 
rate of star formation by a factor of ~ 2, 5, and 15 at z = 2, 3 and 4, respectively. The model is still 
consistent with the global history of light, but overpr edicts the metal mass density at high redshifts 
as sampled by QSO absorbers. In order for such a model to be acceptable, the gas traced by the 
damped Lyman-a systems would have to be physically distinct from the luminous star forming 
regions observed in high-z galaxies, and to be substantially under-enriched in metals compared to 

•^A fitting formula to the stellar birthrate used for Figure 16 is ip(t) = 0.336e" i9 / L6 + 0.0074(1 - 
e -Vo.64) + o.O197t|e _i9 /°- 64 M yr- 1 Mpc- 3 to 1% accuracy. 
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the cosmic mean. 

In the next section we shall compute the expected evolution with cosmic time of the Type la 
and II/Ib,c supernova frequencies for the two star formation histories - for brevity, we shall refer 
to them as the "merging" versus "monolithic collapse" scenarios - discussed above. By focusing on 
the integrated light radiated by the galaxy population as a whole, our approach will not specifically 
address the evolution and the SN rates of particular subclasses of objects, like the oldest ellipticals 
or low-surface brightness galaxies, whose star formation history may have differed significantly from 
the global average. 



2.2. Type la and II(+Ib/c) Supernova Rates 

Single stars with mass > 8M evolve rapidly (< 50 Myr) through all phases of central nuclear 
burning, ending their life as Type II SNe with different characteristics depending on the progenitor 
mass (e.g. Iben & Renzini 1983). For a Salpeter IMF (with lower and upper mass cutoffs of 0.1 
and 125 M Q ), the core-collapse supernova rate can be related to the stellar birthrate according to 



/ 8 125 c?m0(m) _ nnnr7A vy [ if>(t) 



SNRn(t) = ^(t)-355-; ^-^7-^7 = 0.0074 x 



M yr^ 1 Mpc" 



yr" 1 Mpc" 3 . (3) 



/ 1 dmm(j){m) 

It is worth noting at this stage that our model predictions for the frequency of Type II events 
are largely independent of the assumed IMF. This follows from the fact that the rest-frame UV 
continuum emission - which is used as an indicator of the instantaneous star formation rate ip(t) - 
from all but the oldest galaxies is entirely dominated by massive stars on the main sequence, the 
same stars which later give origin to a core collapse SN. 

The specific evolutionary history leading to a Type la event remains instead an unsettled 
question. SN la are believed to result from the explosion of C-0 white dwarfs (WDs) possibly 
triggered by the accretion of material from a companion, the nature of which is still unknown (see 
Ruiz-Lapuente et al. 1997 for a recent review). In a double degenerate (DD) system, for example, 
such elusive companion is another WD: the exploding WD reaches the Chandrasekhar limit and 
carbon ignition occurs at its center (Iben & Tutukov 1984). In the single degenerate (SD) model 
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instead, the companion is a nondegenerate, evolved star that fills its Roche lobe and pours hydrogen 
or helium onto the WD (Whelan & Iben 1973; Iben & Tutukov 1984). While in the latter the clock 
for the explosion is set by the lifetime of the primary star and, e.g., by how long it takes to the 
companion to evolve and fill its Roche lobe, in the former it is controlled by the lifetime of the 
primary star and by the time it takes to shorten the separation of the two WDs as a result of 
gravitational wave radiation. The evolution of the rate depends then, among other things, on 
the unknown mass distribution of the secondary binary components in the SD model, or on the 
distribution of the initial separations of the two WDs in the DD model. 

To shed light into the identification issue and, in particular, on the clock-mechanism for the 
explosion of Type la's, we shall adopt here a more empirical approach. We choose to parametrize 
the rate of Type la's in terms of a characteristic explosion timescale, r - which defines an explosion 
probability per WD assumed to be independent of time - and an explosion efficiency, rj. The former 
accounts for the time it takes in the various models to go from a newly born (primary) WD to the 
SN explosion itself: a spread of "delay" times results from the combination of a variety of initial 
conditions, such as the mass ratio of the binary system, the distribution of initial separations, 
the influence of metallicity on the mass transfer rate and accretion efficiency, etc., folded with 
the various possible mechanisms that lead to a Type la event. The latter simply accounts for the 
fraction of stars in binary systems that, because of unfavorable initial conditions, will never undergo 
a SN la explosion. We consider as possible progenitors all systems in which the primary star has 
an initial mass higher than m m - m = 3M Q (final mass > 0.72 M , Weidemann 1987) and lower 
than m max = 8M (cf Greggio & Renzini 1983): stars less massive than 3M will not produce 
a catastrophic event even if the companion has comparable mass, while stars more massive than 
8 M & will undergo core collapse, generating a Type II explosion. 

With these assumptions the rate of Type la events at any one time will be given by the sum of 
the explosions of all the binary WDs produced in the past that have not had the time to explode 
yet, i.e. 

j] J* ^(t')dt' exp(-^=^U(m)dm 

SNRi a (£) = ' Jme — -P- ^-J , 4 

r J <p(m)dm 
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where m c = max[m m ; n , m(l/)], m(t') = (10 Gyr/i') - 4 is the minimum mass of a star that reaches 
the WD phase at time t', and t m = 10Gyr/m 2 5 is the standard lifetime of a star of mass m (all 
stellar masses are expressed in solar units). For a fixed initial mass m, the frequency of Type 
la events peaks at an epoch that reflect an "effective" delay At = r + t m from stellar birth. A 
prompter (smaller r) explosion results in a higher SN la rate at early epochs. 

3. Comparison with the Data 

Observed rates of SNe are normally given in units of SNu, one SNu corresponding to 1 SN per 
100 years per 10 10 Lbq. Since a galaxy luminosity depends on Hq, there is a factor h 2 involved in 
the inferred SN frequencies. The typical uncertainties on the present-day rates are summarized in 
Table 1, where we have reported the determinations of Cappellaro et al. (1997), Tammann, Lofner, 
& Schroder (1994), and Evans, van den Bergh, & McClure (1989). These agree to within a factor 
of 3 for Type la in early type galaxies and Type Ib/c in spirals, and to within a factor of 2 for Type 
II in late spirals. The "internal" error in each entry is of order 40%. When weighted according to 
the local blue luminosity function by spectral type of Heyl et al. (1997) (which assigns about 28% 
of the total local blue luminosity density to elliptical galaxies, 32% to early spirals, and 40% to 
late spirals), the tabulated frequencies yield a mean - averaged over the entire galaxy population 
- SN la rate of (0.12,0.19,0.12) h 2 50 SNu, and a SN II(+Ib/c) rate of (0.32,0.62,0.60) h 2 m SNu (the 
three values corresponding to the measurements of Cappellaro et al. 1997, Tammann et al. 1994, 
and Evans et al. 1989, respectively). A simple geometric mean of these determinations gives the 
values, SNRi a = 0.14 ±0.06 ft| SNu and SNR n = 0.49 ±0.2 h 2 50 SNu, we shall adopt for comparison 
with the theoretical rates. At higher redshifts, searches for SNe have been pioneered by Hansen et 
al. (1989) and Norgaard-Nielson et al. (1989). Recently, systematic searches of distant SNe (Kim 
et al. 1997; Garnavich et al. 1998; Perlmutter et al. 1998) have provided the first measurement of 
the rate of Type la at z = 0.4, SNRi a = 0.21^;^ h\ SNu (Pain et al. 1997). 

In Figure 2 we show the predicted Type la and II(±Ib/c) rest-frame frequencies as a function 
of redshift. Expressed in SNu, the Type II rate is basically proportional to the ratio between the 
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UV and blue galaxy luminosity densities, and is therefore independent of cosmology. Unlike the SN 
frequency per unit volume, which will trace the evolution of the stellar birthrate, the frequency of 
Type II events per unit blue luminosity is a monotonic increasing function of redshift, and depends 
only weakly on the assumed star formation history. The Type la rates plotted in the figure assume 
characteristic "delay" timescales after the collapse of the primary star to a WD equal to r = 0.3, 1 
and 3 Gyr, which virtually encompass all relevant possibilities. The SN la explosion efficiency was 
left as an adjustable parameter to reproduce the observed ratio of SN II to SN la explosion rates 
in the local universe, SNRn/SNRi a ~ 3.5, 5% < rj < 10% for the adopted models. 

It appears that observational determinations of the SN la rate at z ~ 1 can unambiguously 
identify the appropriate delay time. In particular, we estimate that measuring the frequency of 
Type la events at both z ~ 0.5 and z ~ 1 with an error of 20% or lower would allow one to determine 
this timescale to within about 30%. This kind of observations are by no means prohibitive, and 
these goals could be achieved within a couple of years. In fact, ongoing searches for high-z SNe 
(Perlmutter et al. 1998; Garnavich et al. 1998) are currently able to discover and study about 
a dozen new events per observing session in the redshift range 0.4-1.0, and the observations are 
carried out at a rate of about four sessions a year. Since determining the frequency of SN la with a 
20% uncertainty requires statistics on more than 25 objects per redshft bin, it is clear that, barring 
systematic biases, those rates will soon be know with high accuracy. Also note how, relative to 
the merging scenario, the monolithic collapse model predicts Type la rates (in SNU) that are, in 
the r = 0.3 Gyr case, a factor of 1.6 and 4.9 higher at z = 2 and 4, respectively, with even larger 
factors found in the case of longer delays. Therefore, once such timescale is calibrated through the 
observed ratio SNRi a (z = 0)/SNRi a (z = 1), one should be able to constrain the star formation 
history of the early universe by comparing the predicted SN la rate at z > 2 with the observations. 



4. Conclusions 



We have investigated the link between SN statistics and galaxy evolution. Using recent de- 
terminations of the star formation history of field galaxies from the present epoch to high-z, and 
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a simple model for the evolutionary history of the binary system leading to a Type la event, we 
have computed the theoretical Type la and Type II SN rates as a function of cosmic time. While 
significant uncertainties still remain in these estimates, we believe the calculations presented in 
this Letter offer a first, realistic glimpse to the evolution of the cosmic supernova rates with cosmic 
time. Our main results can be summarized as follows. 

• At the present epoch, the predicted Type II(+Ib/c) frequency appears to match remarkably 
well the observed local value. The obvious caveat is the well known fact that the SN II 
rate is a sensitive function of the lower mass cutoff of the progenitors, mi. Values as low as 
mi = 6 M (Chiosi, Bertelli, & Bressan 1992) or as high as mi = 11 M (Nomoto 1984) have 
been proposed in the literature: adopting a lower mass limit of 6 or 11 M Q would increase 
or reduce our Type II rates by a factor 1.5, respectively. Note also that rates obtained from 
traditional distant (beyond 4 Mpc) sample might need to be increased by a factor of 1.5-2 
because of severe selection effects against Type IPs fainter than My = — 16 (Woltjer 1997). 

• In the interval < z < 1 , the predicted rate of SN la is a sensitive function of the characteristic 
delay timescale between the collapse of the primary star to a WD and the SN event. Accurate 
measurements of SN rates in this redshift range will improve our understanding of the nature 
of SN la progenitors and the physics of the explosions. Ongoing searches and studies of 
distant SNe should soon provide these rates, allowing a universal calibration of the Type la 
phenomenon. 

• While Type la rates at 1 < z < 2 will offer valuable information on the star formation history 
of the universe at earlier epoch, the full picture will only be obtained with statistics on Type 
la and II SNe at redshifts 2 < z < 4 or higher. At these epochs, the detection of Type 
II events must await the Next Generation Space Telescope (NGST). A SN II has a typical 
peak magnitude Mb ~ —17 (e.g. Patat et al. 1994): placed at z = 3, such an explosion 
would give rise to an observed flux of 15 nJy (assuming a flat cosmology with qo = 0.5 and 
Hq = 50/i5o kms -1 Mpc -1 ) at 1.8 /im. At this wavelength, the imaging sensitivity of an 8m 



-12- 



NGST is 1 nJy (10 4 s exposure and 10<7 detection threshold), while the moderate resolution 
(A/ A A = 1000) spectroscopic limit is about 50 times higher (10 5 s exposure per resolution 
element and 10<7 detection threshold) (Stockman et al. 1998). The several weeks period of 
peak rest-frame blue luminosity would be stretched by a factor of (1 + z) to few months. 
Figure 3 shows the cumulative number of Type II events expected per year per 4' x 4' field. 
Depending on the history of star formation at high redshifts, the NGST should detect between 
7 (in the merging model) and 15 (in the monolithic collapse scenario) Type II SNe per field 
per year in the interval 2 < z < 4. The possibility of detecting Type II SNe at z > 5 from 
an early population of galaxies has been investigated by Miralda-Escude & Rees (1997). By 
assuming these are responsible for the generation of all the metals observed in the Lyman-a 
forest at high redshifts, a high baryon density (O&ftJ = 0.1), and an average metallicity of 
0.01Z©, Miralda-Escude & Rees estimate the NGST should observe about 16 SN II per field 
per year with z > 5. Note, however, that a metallicity smaller by a factor ~ 10 compared 
to the value adopted by these authors has been recently derived by Songaila (1997). For 
comparison, the models discussed in this Letter predict between 1 and 10 Type II SNe per 
field per year with z > 4. 

PM acknowledges support from NASA through ATP grant NAG5-4236. MDV thanks the 
ospitality of the STScI, where part of this work was done. 
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Table 1. Present-epoch SN Rates ab 



Hubble Type 


la 


Ib/c 


II 


E-SO 


0.07, 0.25, 0.08 






SOa-Sa 


0.11, 0.12, 0.15 


0.07,0.01, ••• 


0.07, 0.04, • • • 


Sab-Sb 


0.08, 0.12, 0.08 


0.04, 0.07, 0.15 


0.24, 0.34, 0.45 


Sbc-Sd 


0.11, 0.12, 0.05 


0.07, 0.19, 0.18 


0.39, 0.98, 0.58 



a Values listed are in units of /i| SNu (1 SNu= 1 SN per century per 10 10 Lb©)- 
b The first column in each entry shows the value determined by Cappellaro et al. (1997), 
second column by Tammann et al. (1994), third column by Evans et al. (1989). 
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Fig. 1. — Evolution of the observed comoving luminosity density at rest-frame wavelengths of 0.15 
(dotted line), 0.28 (solid line), 0.44 (short-dashed line), 1.0 (long-dashed line), and 2.2 (dot-dashed 
line) /mi. The data points with error bars are taken from Lilly et al. (1996) (filled dots), Connolly 
et al. (1997) (empty squares), Madau et al. (1996) and Madau (1997) (filled squares), Ellis et al. 
(1996) (empty triangles), Treyer et al. (1998) (filled triangle), and Gardner et al. (1997) (empty 
dot). A fiat cosmology with qo = 0.5 and Hq = 50/i5okms _1 Mpc -1 was adopted. The inset in the 
upper-right corner of the plot show the SFR density ( M yr _1 Mpc -3 ) vs. redshift which was used 
as input to the population synthesis code, (a) The model assumes a Salpeter IMF, SMC-type dust 
in a foreground screen (with the scattering term removed), and a universal -Eb-v = 0.1. (b) Test 
case with a much larger star formation density at high redshift than indicated by the HDF dropout 
analysis. The model - designed to mimick a "monolithic collapse" scenario - assumes a Salpeter 
IMF and a dust opacity which increases rapidly with redshift, -Eb-v = 0.011(1 + z) 2 - 2 . 
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Fig. 2. — Predicted Type la and II(+Ib/c) rest-frame frequencies as a function of redshift. The 
rates are normalized to the emitted blue luminosity density. Solid line: SN II rate. Dashed-dotted 
line: SN la rate with r = 0.3 Gyr. Dotted line: SN la rate with r = 1 Gyr. Dashed line: SN la 
rate with r = 3 Gyr. The data points with error bars have been derived from the measurements of 
Cappellaro et al. (1997), Tammann et al. (1994), Evans et al. (1989), and Pain et al. (1997). (a) 
Model predictions for the "merging" scenario of Figure la. (b) Same for the "monolithic collapse" 
scenario of Figure 1 b. 
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Fig. 3. — Predicted cumulative number of Type la and II(+Ib/c) SNe above a given redshift z in 
a 4' x 4' field. Solid line: Type II's. Dashed- dotted line: Type la's with r = 0.3 Gyr. Dotted line: 
Type la's with r = 1 Gyr. Dashed line: Type la's with r = 3 Gyr. The effect of dust extinction 
on the detectability of SNe is negligible in these models, (a) Model predictions for the "merging" 
scenario of Figure la. (b) Same for the "monolithic collapse" scenario of Figure 1 b. 



